The Phenix-AlphaFold webservice: Enabling AlphaFold predictions for use in Phenix.

Advances in machine learning have enabled sufficiently accurate predictions of protein structure to be used in macromolecular structure determination with crystallography and cryo-electron microscopy data. The Phenix software suite has AlphaFold predictions integrated into an automated pipeline that can start with an amino acid sequence and data, and automatically perform model-building and refinement to return a protein model fitted into the data. Due to the steep technical requirements of running AlphaFold efficiently, we have implemented a Phenix-AlphaFold webservice that enables all Phenix users to run AlphaFold predictions remotely from the Phenix GUI starting with the official 1.21 release. This webservice will be improved based on how it is used by the research community and the future research directions for Phenix.

An Engineered Laccase from Fomitiporia mediterranea Accelerates Lignocellulose Degradation.

Laccases from white-rot fungi catalyze lignin depolymerization, a critical first step to upgrading lignin to valuable biodiesel fuels and chemicals. In this study, a wildtype laccase from the basidiomycete Fomitiporia mediterranea (Fom_lac) and a variant engineered to have a carbohydrate-binding module (Fom_CBM) were studied for their ability to catalyze cleavage of ß-O-4' ether and C-C bonds in phenolic and non-phenolic lignin dimers using a nanostructure-initiator mass spectrometry-based assay. Fom_lac and Fom_CBM catalyze ß-O-4' ether and C-C bond breaking, with higher activity under acidic conditions (pH < 6). The potential of Fom_lac and Fom_CBM to enhance saccharification yields from untreated and ionic liquid pretreated pine was also investigated. Adding Fom_CBM to mixtures of cellulases and hemicellulases improved sugar yields by 140% on untreated pine and 32% on cholinium lysinate pretreated pine when compared to the inclusion of Fom_lac to the same mixtures. Adding either Fom_lac or Fom_CBM to mixtures of cellulases and hemicellulases effectively accelerates enzymatic hydrolysis, demonstrating its potential applications for lignocellulose valorization. We postulate that additional increases in sugar yields for the Fom_CBM enzyme mixtures were due to Fom_CBM being brought more proximal to lignin through binding to either cellulose or lignin itself.

Universal parameters of bulk-solvent masks.

The bulk solvent is a major component of biomacromolecular crystals that contributes significantly to the observed diffraction intensities. Accurate modelling of the bulk solvent has been recognized as important for many crystallographic calculations. Owing to its simplicity and modelling power, the flat (mask-based) bulk-solvent model is used by most modern crystallographic software packages to account for disordered solvent. In this model, the bulk-solvent contribution is defined by a binary mask and a scale (scattering) function. The mask is calculated on a regular grid using the atomic model coordinates and their chemical types. The grid step and two radii, solvent and shrinkage, are the three parameters that govern the mask calculation. They are highly correlated and their choice is a compromise between the computer time needed to calculate the mask and the accuracy of the mask. It is demonstrated here that this choice can be optimized using a unique value of 0.6 Šfor the grid step irrespective of the data resolution, and the radii values adjusted correspondingly. The improved values were tested on a large sample of Protein Data Bank entries derived from X-ray diffraction data and are now used in the computational crystallography toolbox (CCTBX) and in Phenix as the default choice.

Accounting for nonuniformity of bulk-solvent: A mosaic model.

A flat mask-based model is almost universally used in macromolecular crystallography to account for disordered (bulk) solvent. This model assumes any voxel of the crystal unit cell that is not occupied by the atomic model is occupied by the solvent. The properties of this solvent are assumed to be exactly the same across the whole volume of the unit cell. While this is a reasonable approximation in practice, there are a number of scenarios where this model becomes suboptimal. In this work, we enumerate several of these scenarios and describe a new generalized approach to modeling the bulk-solvent which we refer to as mosaic bulk-solvent model. The mosaic bulk-solvent model allows nonuniform features of the solvent in the crystal to be accounted for in a computationally efficient way. It is implemented in the computational crystallography toolbox and the Phenix software.

Genome-scale and pathway engineering for the sustainable aviation fuel precursor isoprenol production in Pseudomonas putida.

Sustainable aviation fuel (SAF) will significantly impact global warming in the aviation sector, and important SAF targets are emerging. Isoprenol is a precursor for a promising SAF compound DMCO (1,4-dimethylcyclooctane) and has been produced in several engineered microorganisms. Recently, Pseudomonas putida has gained interest as a future host for isoprenol bioproduction as it can utilize carbon sources from inexpensive plant biomass. Here, we engineer metabolically versatile host P. putida for isoprenol production. We employ two computational modeling approaches (Bilevel optimization and Constrained Minimal Cut Sets) to predict gene knockout targets and optimize the "IPP-bypass" pathway in P. putida to maximize isoprenol production. Altogether, the highest isoprenol production titer from P. putida was achieved at 3.5 g/L under fed-batch conditions. This combination of computational modeling and strain engineering on P. putida for an advanced biofuels production has vital significance in enabling a bioproduction process that can use renewable carbon streams.

Community recommendations on cryoEM data archiving and validation.

In January 2020, a workshop was held at EMBL-EBI (Hinxton, UK) to discuss data requirements for the deposition and validation of cryoEM structures, with a focus on single-particle analysis. The meeting was attended by 47 experts in data processing, model building and refinement, validation, and archiving of such structures. This report describes the workshop's motivation and history, the topics discussed, and the resulting consensus recommendations. Some challenges for future methods-development efforts in this area are also highlighted, as is the implementation to date of some of the recommendations.

Environmental stress mediates groundwater microbial community assembly.

Community assembly describes how different ecological processes shape microbial community composition and structure. How environmental factors impact community assembly remains elusive. Here we sampled microbial communities and >200 biogeochemical variables in groundwater at the Oak Ridge Field Research Center, a former nuclear waste disposal site, and developed a theoretical framework to conceptualize the relationships between community assembly processes and environmental stresses. We found that stochastic assembly processes were critical (>60% on average) in shaping community structure, but their relative importance decreased as stress increased. Dispersal limitation and 'drift' related to random birth and death had negative correlations with stresses, whereas the selection processes leading to dissimilar communities increased with stresses, primarily related to pH, cobalt and molybdenum. Assembly mechanisms also varied greatly among different phylogenetic groups. Our findings highlight the importance of microbial dispersal limitation and environmental heterogeneity in ecosystem restoration and management.

Community recommendations on cryoEM data archiving and validation: Outcomes of a wwPDB/EMDB workshop on cryoEM data management, deposition and validation.

In January 2020, a workshop was held at EMBL-EBI (Hinxton, UK) to discuss data requirements for deposition and validation of cryoEM structures, with a focus on single-particle analysis. The meeting was attended by 47 experts in data processing, model building and refinement, validation, and archiving of such structures. This report describes the workshop's motivation and history, the topics discussed, and consensus recommendations resulting from the workshop. Some challenges for future methods-development efforts in this area are also highlighted, as is the implementation to date of some of the recommendations.

AlphaFold predictions are valuable hypotheses and accelerate but do not replace experimental structure determination.

Artificial intelligence-based protein structure prediction methods such as AlphaFold have revolutionized structural biology. The accuracies of these predictions vary, however, and they do not take into account ligands, covalent modifications or other environmental factors. Here, we evaluate how well AlphaFold predictions can be expected to describe the structure of a protein by comparing predictions directly with experimental crystallographic maps. In many cases, AlphaFold predictions matched experimental maps remarkably closely. In other cases, even very high-confidence predictions differed from experimental maps on a global scale through distortion and domain orientation, and on a local scale in backbone and side-chain conformation. We suggest considering AlphaFold predictions as exceptionally useful hypotheses. We further suggest that it is important to consider the confidence in prediction when interpreting AlphaFold predictions and to carry out experimental structure determination to verify structural details, particularly those that involve interactions not included in the prediction.

Improved joint X-ray and neutron refinement procedure in Phenix.

Neutron diffraction is one of the three crystallographic techniques (X-ray, neutron and electron diffraction) used to determine the atomic structures of molecules. Its particular strengths derive from the fact that H (and D) atoms are strong neutron scatterers, meaning that their positions, and thus protonation states, can be derived from crystallographic maps. However, because of technical limitations and experimental obstacles, the quality of neutron diffraction data is typically much poorer (completeness, resolution and signal to noise) than that of X-ray diffraction data for the same sample. Further, refinement is more complex as it usually requires additional parameters to describe the H (and D) atoms. The increase in the number of parameters may be mitigated by using the `riding hydrogen' refinement strategy, in which the positions of H atoms without a rotational degree of freedom are inferred from their neighboring heavy atoms. However, this does not address the issues related to poor data quality. Therefore, neutron structure determination often relies on the presence of an X-ray data set for joint X-ray and neutron (XN) refinement. In this approach, the X-ray data serve to compensate for the deficiencies of the neutron diffraction data by refining one model simultaneously against the X-ray and neutron data sets. To be applicable, it is assumed that both data sets are highly isomorphous, and preferably collected from the same crystals and at the same temperature. However, the approach has a number of limitations that are discussed in this work by comparing four separately re-refined neutron models. To address the limitations, a new method for joint XN refinement is introduced that optimizes two different models against the different data sets. This approach is tested using neutron models and data deposited in the Protein Data Bank. The efficacy of refining models with H atoms as riding or as individual atoms is also investigated.

Large bi-axial tensile strain effect in epitaxial BiFeO3 film grown on single crystal PrScO3.

A BiFeO3 film is grown epitaxially on a PrScO3 single crystal substrate which imparts ~ 1.45% of biaxial tensile strain to BiFeO3 resulting from lattice misfit. The biaxial tensile strain effect on BiFeO3 is investigated in terms of crystal structure, Poisson ratio, and ferroelectric domain structure. Lattice resolution scanning transmission electron microscopy, precession electron diffraction, and X-ray diffraction results clearly show that in-plane interplanar distance of BiFeO3 is the same as that of PrScO3 with no sign of misfit dislocations, indicating that the biaxial tensile strain caused by lattice mismatch between BiFeO3 and PrScO3 are stored as elastic energy within BiFeO3 film. Nano-beam electron diffraction patterns compared with structure factor calculation found that the BiFeO3 maintains rhombohedral symmetry, i.e., space group of R3c. The pattern analysis also revealed two crystallographically distinguishable domains. Their relations with ferroelectric domain structures in terms of size and spontaneous polarization orientations within the domains are further understood using four-dimensional scanning transmission electron microscopy technique.

Simultaneous carbon catabolite repression governs sugar and aromatic co-utilization in Pseudomonas putida M2.

Pseudomonas putida have emerged as promising biocatalysts for the conversion of sugars and aromatic compounds obtained from lignocellulosic biomass. Understanding the role of carbon catabolite repression (CCR) in these strains is critical to optimize biomass conversion to fuels and chemicals. The CCR functioning in P. putida M2, a strain capable of consuming both hexose and pentose sugars as well as aromatic compounds, was investigated by cultivation experiments, proteomics, and CRISPRi-based gene repression. Strain M2 co-utilized sugars and aromatic compounds simultaneously; however, during cultivation with glucose and aromatic compounds (p-coumarate and ferulate) mixture, intermediates (4-hydroxybenzoate and vanillate) accumulated, and substrate consumption was incomplete. In contrast, xylose-aromatic consumption resulted in transient intermediate accumulation and complete aromatic consumption, while xylose was incompletely consumed. Proteomics analysis revealed that glucose exerted stronger repression than xylose on the aromatic catabolic proteins. Key glucose (Eda) and xylose (XylX) catabolic proteins were also identified at lower abundance during cultivation with aromatic compounds implying simultaneous catabolite repression by sugars and aromatic compounds. Reduction of crc expression via CRISPRi led to faster growth and glucose and p-coumarate uptake in the CRISPRi strains compared to the control, while no difference was observed on xylose+p-coumarate. The increased abundances of Eda and amino acid biosynthesis proteins in the CRISPRi strain further supported these observations. Lastly, small RNAs (sRNAs) sequencing results showed that CrcY and CrcZ homologues levels in M2, previously identified in P. putida strains, were lower under strong CCR (glucose+p-coumarate) condition compared to when repression was absent (p-coumarate or glucose only).IMPORTANCEA newly isolated Pseudomonas putida strain, P. putida M2, can utilize both hexose and pentose sugars as well as aromatic compounds making it a promising host for the valorization of lignocellulosic biomass. Pseudomonads have developed a regulatory strategy, carbon catabolite repression, to control the assimilation of carbon sources in the environment. Carbon catabolite repression may impede the simultaneous and complete metabolism of sugars and aromatic compounds present in lignocellulosic biomass and hinder the development of an efficient industrial biocatalyst. This study provides insight into the cellular physiology and proteome during mixed-substrate utilization in P. putida M2. The phenotypic and proteomics results demonstrated simultaneous catabolite repression in the sugar-aromatic mixtures, while the CRISPRi and sRNA sequencing demonstrated the potential role of the crc gene and small RNAs in carbon catabolite repression.

A combinatorial droplet microfluidic device integrated with mass spectrometry for enzyme screening.

Mass spectrometry (MS) enables detection of different chemical species with a very high specificity; however, it can be limited by its throughput. Integrating MS with microfluidics has a tremendous potential to improve throughput and accelerate biochemical research. In this work, we introduce Drop-NIMS, a combination of a passive droplet loading microfluidic device and a matrix-free MS laser desorption ionization technique called nanostructure-initiator mass spectrometry (NIMS). This platform combines different droplets at random to generate a combinatorial library of enzymatic reactions that are deposited directly on the NIMS surface without requiring additional sample handling. The enzyme reaction products are then detected with MS. Drop-NIMS was used to rapidly screen enzymatic reactions containing low (on the order of nL) volumes of glycoside reactants and glycoside hydrolase enzymes per reaction. MS "barcodes" (small compounds with unique masses) were added to the droplets to identify different combinations of substrates and enzymes created by the device. We assigned xylanase activities to several putative glycoside hydrolases, making them relevant to food and biofuel industrial applications. Overall, Drop-NIMS is simple to fabricate, assemble, and operate and it has potential to be used with many other small molecule metabolites.

Quantitative Analysis of The High-Yield Hydrolysis of Kelp by Laminarinase and Alginate Lyase.

Kelp is an abundant, farmable biomass-containing laminarin and alginate as major polysaccharides, providing an excellent model substrate to study their deconstruction by simple enzyme mixtures. Our previous study showed strong reactivity of the glycoside hydrolase family 55 during hydrolysis of purified laminarin, raising the question of its reactivity with intact kelp. In this study, we determined that a combination of a single glycoside hydrolase family 55 ß-1,3-exoglucanase with a broad-specificity alginate lyase from the polysaccharide lyase family 18 gives efficient hydrolysis of untreated kelp to a mixture of simple sugars, that is, glucose, gentiobiose, mannitol-end glucose, and mannuronic and guluronic acids and their soluble oligomers. Quantitative assignments from nanostructure initiator mass spectrometry (NIMS) and 2D HSQC NMR spectroscopy and analysis of the reaction time-course are provided. The data suggest that binary combinations of enzymes targeted to the unique polysaccharide composition of marine biomass are sufficient to deconstruct kelp into soluble sugars for microbial fermentation.

Alkaline-SDS cell lysis of microbes with acetone protein precipitation for proteomic sample preparation in 96-well plate format.

Plate-based proteomic sample preparation offers a solution to the large sample throughput demands in the biotechnology field where hundreds or thousands of engineered microbes are constructed for testing is routine. Meanwhile, sample preparation methods that work efficiently on broader microbial groups are desirable for new applications of proteomics in other fields, such as microbial communities. Here, we detail a step-by-step protocol that consists of cell lysis in an alkaline chemical buffer (NaOH/SDS) followed by protein precipitation with high-ionic strength acetone in 96-well format. The protocol works for a broad range of microbes (e.g., Gram-negative bacteria, Gram-positive bacteria, non-filamentous fungi) and the resulting proteins are ready for tryptic digestion for bottom-up quantitative proteomic analysis without the need for desalting column cleanup. The yield of protein using this protocol increases linearly with respect to the amount of starting biomass from 0.5-2.0 OD*mL of cells. By using a bench-top automated liquid dispenser, a cost-effective and environmentally-friendly option to eliminating pipette tips and reducing reagent waste, the protocol takes approximately 30 minutes to extract protein from 96 samples. Tests on mock mixtures showed expected results that the biomass composition structure is in close agreement with the experimental design. Lastly, we applied the protocol for the composition analysis of a synthetic community of environmental isolates grown on two different media. This protocol has been developed to facilitate rapid, low-variance sample preparation of hundreds of samples and allow flexibility for future protocol development.

Overall protein structure quality assessment using hydrogen-bonding parameters.

Atomic model refinement at low resolution is often a challenging task. This is mostly because the experimental data are not sufficiently detailed to be described by atomic models. To make refinement practical and ensure that a refined atomic model is geometrically meaningful, additional information needs to be used such as restraints on Ramachandran plot distributions or residue side-chain rotameric states. However, using Ramachandran plots or rotameric states as refinement targets diminishes the validating power of these tools. Therefore, finding additional model-validation criteria that are not used or are difficult to use as refinement goals is desirable. Hydrogen bonds are one of the important noncovalent interactions that shape and maintain protein structure. These interactions can be characterized by a specific geometry of hydrogen donor and acceptor atoms. Systematic analysis of these geometries performed for quality-filtered high-resolution models of proteins from the Protein Data Bank shows that they have a distinct and a conserved distribution. Here, it is demonstrated how this information can be used for atomic model validation.

Efficient structure-factor modeling for crystals with multiple components.

Diffraction intensities from a crystallographic experiment include contributions from the entire unit cell of the crystal: the macromolecule, the solvent around it and eventually other compounds. These contributions cannot typically be well described by an atomic model alone, i.e. using point scatterers. Indeed, entities such as disordered (bulk) solvent, semi-ordered solvent (e.g. lipid belts in membrane proteins, ligands, ion channels) and disordered polymer loops require other types of modeling than a collection of individual atoms. This results in the model structure factors containing multiple contributions. Most macromolecular applications assume two-component structure factors: one component arising from the atomic model and the second one describing the bulk solvent. A more accurate and detailed modeling of the disordered regions of the crystal will naturally require more than two components in the structure factors, which presents algorithmic and computational challenges. Here an efficient solution of this problem is proposed. All algorithms described in this work have been implemented in the computational crystallography toolbox (CCTBX) and are also available within Phenix software. These algorithms are rather general and do not use any assumptions about molecule type or size nor about those of its components.